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determine the mass of the top quark to be A4top = 161 ± 17(stat.) ± 10(syst.) GeV/c^. In addition 
we measure a ti production cross section of pb (where Mtop = 175 GeY/c? has been assumed 

for the acceptance estimate). 

PACS numbers: 14.65.Ha, 13.85.Qk, 13.85.Ni 
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We report here on a measurement of the tt produc- 
tion cross section and top quark mass in the dilep- 
ton channel with the CoUider Detector at Fermilab 
(CDF). The data sample corresponds to a total inte- 
grated luminosity of 109 ± 7 pb~^. This analysis con- 
siders dilepton events originating predominantly from 
tt W+bW-b {(.+vb){(.-vb), with £ = e or A 
subset of these events corresponding to a data sample of 
67 pb^^ supported the discovery of the top quark 

Recently both the CDF and DO collaborations have 
updated their measurements of the ti production cross 
section and top quark mass using the "lepton 
plus jets" channel: tt W+bW~b -> {i+vb){qq'b), with 
(. = e ov fi, which has larger statistics. The consistency of 
the measurements presented here with those in the lepton 
plus jets channel is an important confirmation that these 
two orthogonal sets of events both originate from the 
same heavy top production process. 

A description of the CDF detector can be found in Ref. 
[0. The coordinate system and various quantities used 
throughout this paper are defined in The momenta 
of the charged leptons are measured with the central 
tracking chamber in a 1.4-T superconducting solenoidal 
magnet. Electromagnetic and hadronic calorimeters sur- 
rounding the tracking chambers are used to identify 
and measure the energies of electrons and jets. Muons 
are identified with drift chambers located outside the 
calorimeters. A three-level trigger selects high transverse 
momentum (Pt) electrons and muons. 

The event selection is very similar to the previous anal- 
ysis described in Ref. ||^,^. We require two high-PT- 
(Pt > 20 GeV/c), oppositely charged leptons (e or /i) 
in the central pseudorapidity region {\ri\ < 1.0) with at 
least one of them well isolated from nearby tracks and 
calorimeter activity. We reject Z — > £~^£~X events by 
requiring the dilepton invariant mass, Mee or Af^^, to 
be outside the interval 75-105 GeV/c^. If there is a 
high transverse energy {Et > 10 GeV) isolated photon 
present, the event is removed if consistent with a radia- 
tive Z decay. We also require that there be at least two 
jets with measured Ex > 10 GeV in the pseudorapidity 
range < 2.0, as expected from the presence of two 
b quarks in a tt event. The signature of the two neu- 
trinos in the decay is missing energy transverse to the 
beam direction (^t)', we require |^t| > 25 GeV. The 
is corrected for non-uniformities in calorimeter response 
and absolute energy scale, and for high-Pr muons. To 
ensure that the is not due to mismeasurements of the 
energies of the leptons or jets, we require \^t\ > 50 GeV 
if the azimuthal angle between the direction of the 
vector and the nearest lepton or jet(j), A(/)(^t,^ or j), 
is less than 20° 

The distribution of A(/)(^t, ^ or j) versus Ifrl is shown 
in Figure ^ for dilepton events that pass the invariant 
mass and two-jet requirements. Superimposed is the 
expected distribution from the HERWIG ti Monte 



Carlo program for Mtop = 175 GeV/c^ followed by a de- 
tector simulation|^. We find nine candidate events in the 
signal region: seven e/i, one /x/x and one ee event. 

The identification efficiencies for single leptons are 
measured from Z — > £~^i~ events in the data and are 
found to be 91% for muons and 83% for electrons [ p^ . 
The acceptance for tt decays to pass all the selection cri- 
teria, including lepton identification, is the average of 
the results from the HERWIG and PYTHIA (l|] Monte 
Carlo programs followed by a detector simulation. We 
find for Mtop = 175 GeV/c^ that (0.74 ± 0.08)% of all 
tt decays pass the above dilepton selection criteria. The 
uncertainty is dominated by the differences between the 
event generators and by the systematic uncertainties in 
the detector simulation. The acceptance increases by 
35% as Mtop increases from 150 to 200 GeV/c^. In 
(86 ± 2)% of the dilepton events passing all selection cri- 
teria both leptons come directly from the decays of the 
W bosons. The remainder consists mostly of events in 
which one of the W bosons decays to a t lepton, which 
in turn decays to an electron or muon. Of the dilepton 
events, we expect (58 ± 2)% to be e/i, (27 ± 1)% fifi and 
(15 ± 1)% ee, where the uncertainty is statistical only. 
Using a tt production cross section value of 5.5 pb, con- 
sistent with recent theoretical calculations [|l^ that as- 
sume Mtop=l75 GeV/c^, we expect to observe 4.4 signal 
events. 

As backgrounds we consider standard model processes, 
other than tt, that can result in dilepton final states. The 
main sources are Drell-Yan {Z*/"f ee,fifi), Z tt 
and WW production. If these events contain additional 
jets from QCD radiation plus fix, either from real neu- 
trinos or from energy mismeasurement, they may satisfy 
our selection criteria. Background contributions from ra- 
diative Z bosons and from bb, WZ, ZZ and Wbb pro- 
duction are estimated to be small. Additional sources 
are processes with a real lepton and a jet or a track fak- 
ing a second lepton, and processes in which mismeasured 
muon tracks can result in an overestimate of the fx in the 
event. This latter background is not relevant for electrons 
because the electron energy is measured in the calorime- 
ters. We estimate the background from Drell-Yan pro- 
duction, fake leptons and mismeasured tracks from the 
data; the other backgrounds are calculated using Monte 
Carlo simulations. The background contributions from 
the different sources are listed in Table I. The errors 
include both systematic and statistical uncertainties. Of 
the 2.4 ±0.5 total background events estimated, 0.8 ±0.2 
are expected in the e/i channel, which docs not have a 
contribution from Drell-Yan production, the dominant 



^The top mass value of 175 GeV/c^ is used to be consistent 
with CDF's most precise mass measurement, that from the 
lepton plus jets channel |fn| . 
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TABLE I. Expected dilepton backgrounds in 109 ± 7 pb 



Background type Expected ^ of events 



Urell- Van 


U.Di 


1 

± 


n on 
0.3U 


Z ^ TT 


0.59 


± 


0.14 


Fake leptons 


0.37 


± 


0.23 


WW 


0.36 


± 


0.11 


Mismeasured muon tracks 


0.3 


± 


0.3 


bh 


0.05 


± 


0.03 


Other(radiative Z, Whb, WZ, ZZ) 


0.1 


± 


0.1 


Total 


2.4 


± 


0.5 



background source in the ee and /i/x channels. 

When relaxing the two-jet requirement in the data se- 
lection, we find eight dilepton events with zero jets and 
eleven with one jet, while we expect 8±2 and 7±2 events, 
respectively, from both background and tt. The ti contri- 
bution is small; using the measured ti production cross 
section in this channel (see below) we expect 0.03 ± 0.02 
dilepton events in the zero jet sample and 1.1 ±0.5 events 
with one jet. Of the eleven events in the data with one jet, 
four are dimuons, three of which show the expected char- 
acteristics of the mismeasured muon track background. 
The expectation from this background is 1.4± 1.5 events. 
Five of the remaining seven events are e/i,. One of these 
five events has a jet tagged as a 6 quark (see below). The 
expected number of ti events in the one jet sample with 
the jet tagged as a & quark is about 0.2. 

We find two events in the data that satisfy all selection 
criteria, except for the opposite sign requirement on the 
lepton charge. We expect 0.37 ± 0.23 same sign dilep- 
ton plus two-jet events from fake lepton background and 
0.24 ±0.11 from ti. Without any jet requirements these 
estimates are 2.6±1.9 and 0.29±0.13, respectively. There 
are no events in the data with two same sign leptons and 
or 1 jets that satisfy all the other selection criteria. 

Events originating from ti decays are characterized by 
the presence of jets originating from h quarks. Four of the 
nine candidate events have one jet tagged as a 6 quark, 
as evidenced by the presence of a secondary vertex in the 
silicon vertex detector (SVX tag) Two of these four 
jets are also tagged by a soft lepton from the semileptonic 
decay of the b quark (SLT tag) ||^. No jets are found 
to be tagged by the SLT method alone. In a sample of 
seven ti events, 4.3 ± 0.4 jets are expected to be tagged 
by at least one of the two tagging methods. If all nine 
candidates were from background processes, the number 
of expected ^-tagged jets would be 0.7 ± 0.2. 

Using the nine observed dilepton events, 2.4 ± 0.5 of 
which are estimated to be background, an overall accep- 
tance of (0.74±0.08)% (for Mtop = 175 GeV/c^), and an 
integrated luminosity of 109 ± 7pb~^, we calculate the ti 
production cross section to be 8.2t3;4 pb. 

With ti dilepton decays, reconstructing the top quark 



mass from the measured final state is a kinematically 
under-constrained problem due to the presence of two 
neutrinos. 

To increase the purity of the sample for the mass anal- 
yses, we require Ht > 170 GeV. Ht is defined as the sum 
of the Pt's of the two leptons [Et for electrons), the Et's, 
of the two highest Et jets and the |^t|- For the nine can- 
didate events the Ht distribution is shown in Figure |^, 
together with the expectations from a Monte Carlo cal- 
culation for Mtop ~ 175 GeV/c^ and background. The 
efficiency of the Ht > 170 GeV cut is about 95% for 
Mtop = 175 GeV/c^. After this cut, eight of the nine 
dilepton candidate events survive and the background is 
reduced to 1.3 ± 0.3 events. 

We present two methods to determine the top quark 
mass. The first method takes advantage of the corre- 
lation between the energy of the jets from b quarks in 
top quark decays and the top quark mass. From Monte 
Carlo simulation of ti decays we find that the mean en- 
ergy of the two highest Et jets increases linearly with 
Mtop with a slope of 0.5. We obtain the most probable 
value of Mtop for our set of events by comparing the ob- 
served distribution of jet energies with Monte Carlo dis- 
tributions (templates) for different top quark masses. We 
use the two highest Et jets in the event. The jet energies 
are corrected for non-uniformities in calorimeter response 
and absolute energy scale. The templates for ti decays 
are obtained from large samples of Monte Carlo events 
generated with HERWIG, and reconstructed with the 
same analysis programs as the data. Templates are also 
obtained for the background processes discussed above. 
The top quark mass is estimated by performing a max- 
imum likelihood fit of the jet energy distribution from 
the data to a combination of ti and background tem- 
plates. The amount of background is constrained to the 
expected value within its uncertainty. From the resulting 
likelihood values. Cm, for each assumed top quark mass 
TO, the negative logarithms, — lii{Cm), are fit with a third 
order polynomial and the value of the mass correspond- 
ing to the minimum of the polynomial is obtained. The 
statistical uncertainties are determined by the values of 
the mass that give an increase of 0.5 in — ln[Cm) from 
the minimum. 

We have tested the fitting method with many Monte 
Carlo samples, with Mtop in the range between 100 and 
240 GeV/c^, and with the same statistics as the data. In- 
cluded in these samples is the expected number of back- 
ground events. The mean of the distribution of the es- 
timated top quark mass for each sample agrees with the 
generated mass within 2 GeV/c^. The average statis- 
tical uncertainties derived from the likelihood fits are 
consistent with the spread of the distributions, about 
21 GeV/c^, and are not strongly dependent on the top 
quark mass. 

Figure ||(a) shows the jet energy distribution of the 
two highest Et jets for the eight events in the data over- 
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layed with a template obtained from a combination of 
tt Monte Carlo {Mtop = 160 GeV/c^) and background 
in the ratio of 6.7 to 1.3. The distribution of the back- 
ground alone is also shown. The inset shows the poly- 
nomial fit to the — hi{Crn) values versus top quark mass. 
With this method we measure Mtop — 159 ± 23(stat.) ± 
ll(syst.)GeV/c2. 

The second method uses the invariant mass Ma of 
the charged lepton £ and the b quark. In the decay of 
top quarks {t Wb), the energy of the b quark in the 
rest frame of the W boson, Et, is a constant {M^^^ = 
M^ + M^ + 2MwEb). With the subsequent semileptonic 
decay W £v, Eh can be obtained from the invariant 
mass Mill of the charged lepton i and the b quark, and 
their opening angle cos 9 a ■ For a sample of tt dilepton 

events, Mtop is related to Ma as Mfop = ^w' + T^^^^^' 
where terms that include lepton and b quark masses have 
been neglected. (M|j) and {cos da) are the mean values 
of A/^jj and cos 9a in the sample. In the standard model 
tree-level calculation, (cos 9a) = M^/(MjQp-|-2M^), re- 
suiting in Ml^ = (M|) + M^ + iM'^ (M^) + {M'^,)^ . 

In each dilepton tt candidate, one of the two lepton- 
jet combinations {£~^ji,i~j2), {^~^j2,£~ji) corresponds to 
the correct {£^b,£~b) assignment. We select the combi- 
nation with the smallest sum of invariant masses, result- 
ing in two values of A/^^^^j„ per event. The probability of 
picking up the correct combination is in the range of 55% 
to 75%, depending on Mtop- The distribution of 
for the eight candidate events is shown in Figure ^(b) 
together with the expectation for Mtop — 160 GeV/c^. 
The mean value of M^,^.^, (7.2 ± 1.9) x 10^ GeV/c^, is 
used with the linear mapping function {M'^^) — Cq + 
Ci{Ml^J (with Co = (-2.90 ±0.26) x lO^, Ci = 
1.57 ±0.03) to obtain the value of (A^|f,) corresponding to 
standard model ti production and decay. The mapping 
function is determined from ti Monte Carlo events gen- 
erated with HERWIG and a CDF detector simulation. 
It accounts for selection biases, incorrect lepton-jet com- 
binations, and jet energy mismeasurements. We obtain 
Mtop = 163 ± 20(stat.) ± 9(syst.) GeV/c^. 

We combine the two mass results taking into account 
the correlations in the uncertainties to obtain a single 
value for the top quark mass, Mtop — 161 ± 17(stat.) ± 
10(syst.) GeV/c^. The major contributions to the sys- 
tematic uncertainty are due to the uncertainty in the jet 
energy scale and in the shape of the background distri- 
butions. 

In conclusion, in the dilepton channel we find nine can- 
didate events consistent with originating from tt produc- 
tion. The estimated background is 2.4 ± 0.5 events. We 
measure the ti production cross section to be 8.2^3'4 pb 
for Mtop = 175 GeV/c^. The measured top quark mass 
is Mtop = 161 ± 17(stat.) ± 10(syst.) GeV/c^ consistent 
with the CDF measurement in the lepton plus jets chan- 
nel of 175.9 ± 4.8(stat.) ± 4.9(syst.) GeV/c^ Q. 
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FIG. 1. Azimuthal angle between and the nearest lepton or jet, versus \^t\ for events with two leptons and two jets. 
The dashed hne represents the cut. The small dots are for tt Monte Carlo for Mtop = 175 GeV/c^ and correspond to an 
integrated luminosity of about 24 fb~^. The larger symbols represent the data. 
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FIG. 2. Comparison of Ht for the candidate events (histogram) and the expectation from tt production {Mtop = 175 GeV/c^) 
plus background (lighter shaded area). The background distribution alone is represented by the darker shaded area. 
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FIG. 3. (a) Jet energy distribution of the two highest Et jets for the dilepton events (histogram). Superimposed is the same 
distribution for tt Monte Carlo {Mtc,p—160 GeV/c^) plus background (lighter shade), and the background alone (darker shade). 
The tt plus background distribution is normalized to the data. In the inset we show the — ln{Cm) fit as a function of Mtop 
(the minimum has been offset to be at zero), (b) Distribution of Mf^,^-^ (histogram). Superimposed is the expectation from 
ti Monte Carlo {Mtop=160 GeV/c^) plus background (lighter shade), and background alone (darker shade). 
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